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(N : ABSTRACT 

We explore the evolution of the large-scale anisotropy in the velocity field 
caused by the gravitational pancaking effect assuming a ACDM universe. The 
On| ■ Millennium Run halo catalogs at four different redshifts, 2; = 0, 0.5, 1 and 

z = 2 are analyzed to find that the pancaking effect starts to intervene the 



> 



^ . hierarchical structure formation at redshift z = 2 when a characteristic pancake 



X 
J3 



scale is around 3 h ^Mpc. It is also clearly shown how the degree and scale of the 



O , pancaking effect changes with time. An analytic model based on the Zel'dovich 

O . approximation is presented to explain quantitatively the evolution of the velocity- 

pancake alignment. A cosmological implication of our finding and a possibility 
of detecting a signal in real universe are discussed. 

O 

c/3 . Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

The observed filamentary distribution of the large scale structure in the universe, which 
are often called the co smic web, have so far motiva t ed pl enty of works. After the ingenious 



breakthrough made by iBond. Kofman. &: Pogosyanl (119961 ) , a theoretical framework has been 



provided within which the characteristic properties of the cosmic web can be explained 
qualitatively in terms of the tidal infiuences in the universe dominated by the cold dark 
matter (CDM). 

The influence of the tidal forces on the orientations and distributions of the large scale 
structure has been a target of many fruitful studies which either numerically or observation- 
ally proved that the tidal effect is strongest near the pancakes. 



A pancake refers to the two-dimensional sheet-like structure in the universe, which 
form through gravitational collapse along the local tidal field. The first order Lagrangian 
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perturbation theory, the Zel'dovich approximation (jZel'dovichl llQTOl ) . is probably the most 
popular model for the pancake formation which was originally suggested assuming the hot 
dark matter (HDM) dominated universe. The cosmic web theory, however, showed brilliantly 
that the presence of pancakes are also very possible in the CDM-dominated universe due 
to the large scale coherence of the tidal field. N-body simulation s of CDM particles in deed 
demonstrated that the first collapsed objects look like pancakes (jShandarin et al.l 119951 ). 



The large-scale coherence of the tidal field that is responsible for the formation of pan- 
cakes in turn causes the dark matter halos near or on the pancakes to have strong spatial 



correlations in the orientations (lAragon-Calvo et al.l 120071 : iBrunino et al.l 120071 : iHahn et al. 



20071 ). The reason that the degree of the tidally induced alignments is highest near the 
pancakes can be understood given the fact that the pancakes are the first collapsed objects, 
being still in linear regime while the filaments are more severely modified by the nonlinear 
process. 

The formation of pancakes induces not only the spatial correlations in the halo orienta- 



tions but also the anisotropy in the velocity field. In our previous paper (INoh and Led 12006 



hereafter, NL06), we have for the first time detected a 5cr significant signal of the velocity 
anisotropy induced by the pancaking effect in the high-resolution N-body simulation. In 
spite that the signal is found to be very weak because of the non-conserved nature of the 
velocity field, the velocity anisotropy holds a crucial key to understanding the pancaking 
effect on the structure formation. 

In the light of the first detection of NL06, a couple of questions naturally arise: When 
does the pancaking effect influence most dominantly the formation of dark halos?; How does 
the pancake scale change with time?; What are the implications and consequences of the 
pancaking effect on the structure formation?; Is the standard theory capable of quantifying 
the velocity anisotropy?. Our goal here is to provide answers to these questions. 

We organize this paper as follows. In §2, we introduce the analyzing method to in- 
vestigate the evolution of the pancaking effect in the Millennium Simulation. We describe 
analytic modeling and compare it with the results obtained in the simulation catalog in §3. 
Finally, we discuss and conclude the results. 



SIGNALS FROM SIMULATION 



To investigate the evolution of local pancaking effect, we use the Millennium Run simu- 
lation halo catalogs for the concordance ACDM cosmology (jSpringel et alll2005l ). A periodic 
box of this simulation is 500 h~^Mpc on each side and the cosmological parameters are 
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= 0.25, l^A = 1 — fim, h = 0.73, and as = 0.9. We choose four different redshift cata- 
logs, approximately z = 0, 0.5, 1, and 2 among 64 output times of the simulation. Then 
we analyze the alignments of dark-matter halo velocities to the local pancake plane in each 
catalog. Basically, to find a local pancake, we follow the methodology suggested by NL06. 

As a first step in measurement, we select the halos which contain more than 50 particles 
(the left column of Table [1]) since those halos in Millennium Run catalog which have too 
low particle numbers (less than 30) are poorly sampled and severely contaminated by the 
numerical noise. In other words, only those halos which have enough particles to define a 
halo density profile are reliable (V. Springel in private communication). Thus, we set the 
particle number cut-off at 50. Next, among the selected halos, w e pick out field halos w hich 



are still in quasi- li near regime by using the method suggested by lEl-Ad &: PiranI (119971 ) and 



Hoyle fc Vogeleyi (|2002| ). The number of identified field halos is shown in the right column 
of Table [H The field halos are suited for our purpose because they may not be seriously 
affected by gravitational forces of many other halos, retaining their initial conditions at the 
moment that they collapsed. 

Using only those field halos, we determine local pancakes using the practical methods 
proposed by NL06. First, We find two nearest field halos for a field halo and define a 
local pancake plane enclosing all the three halos. Here, the local pancakes are found by 
changing the criterion distance Rc from h~^Mpc to 17 h~^Mpc which is a lower limit 
of Ri, the displacement vector to the first nearest halo. In other words, Rc satisfies the 
following condition: -Rc < -Ri < R2 where R2 is the displacement vector to the second 
nearest neighbor. 

With the determined local pancakes, then we calculate the alignment angles between a 
halo velocity vector and a vector normal to the local pancake plane, cos^. Finally, we obtain 
the probability density distribution of cos 9, p(cos 9), by counting the number of halos in each 
bin of cos 6'. We note that in order to examine the dependence of the particle number cut-off 
which is used when we determine the halos, we recalculate the velocity-pancake alignment 
for those halos which include more than 30 particles though low-particle number halos suffer 
from simulation noise. From this calculation, we obtain almost perfectly consistent results 
with the results of the catalog in which halos include more than 50 particles as can be seen 
in Fig. [1] (particle number cut-off is 50) and Fig. [2] (particle number cut-off is 30). Thus, 
the velocity-pancake alignment is not affected by the choice of particle number cut-off. 

Figure [T] shows p{cos9) with Poissonian error. The three columns correspond to three 
different cases of 2; = 0, 1, and 2, respectively, and the three rows to three different scales of 
Rc = 0, 3, and 6, respectively. The horizontal dotted line in each panel represents a random 
distribution. As can be seen in the right column, when i?c is 0, the halo velocity tends to be 
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slightly aligned with the normal to the local pancake plane at z = 2. This tendency becomes 
stronger as Rc increases to Gh^^Mpc. It can be interpreted that the local pancaking effect 
arises at least around z ~ 2. Meanwhile, at z = 1, for the case of Rc = 0, there is no 
alignment. Then local pancaking effect starts to be visible near R^ = 3/i~^Mpc. Moreover, 
as one can see in the left column, when Rc = at z = 0, the halo velocity does not tend to be 
aligned with the normal to the local pancake plane any more. Rather it tends to lie on the 
plane. The alignment signal begins to be seen when Rc is approximately 6/i~^Mpc. We infer 
from these phenomena that the pancaking effect disappears since the gravitational attraction 
forces from the close neighbor halos affect halo velocities rather than the initial pancaking 
effect. Therefore, we find the velocity-pancake alignments are strongest at z ~ 2, indicating 
that the pancaking effect starts at 2; ^ 2. We refer that we do not make more analysis of the 
velocity-pancake alignments for higher redshift to find a redshift at which pancaking effect 
starts since the number of halos decreases rapidly as redshift increases. Even if the value 
of rj itself might increase, the value of cr,, which is inversely proportional to the number of 
halos would increase. Thus, it would be hard to determine precisely at which redshift both 
the values of 77 and cr^ reach maximums. 

Figure [3] plots the average of cos 6* at four different z = 0, 0.5, 1, and 2 as a function of 
linear characteristic pancake scale, Lp. We refer that Lp is defined as the average distance 
between the center of mass (CM) for three halos on the local pancake plane to all the three 
halos, which was proposed by NL06. The horizontal dotted line in each panel represents 
(cos 6) of random distribution of the angles. As one can see in each panel, there is a specific 
range of Lp where {cos 6) is relatively higher than those at the other pancake scales. This 
specific range implies local pancaking effect significantly occurs on specific pancake scale. 
Comparing to panels each other, overall (cos 9) increases, its slope is steeper and the specific 
range shifts to small Lp when going to high redshift. This tendency can be clearly seen when 
we see that the regions of filled with oblique lines, which represent the area consisting of a 
few of the highest values, move to small Lp as going to high redshift. More quantitatively, 
the peak value is approximately 5.035 x 10~^ at z = while the one is 5.055 x 10"^ at z = 2. 
In addition, the peak position of histogram is approximately 2 h~^Mpc at z = whereas the 
one is 8 h~^Mpc at z = 2. This gradual change with redshifts can be understood the same 
way as Fig. [H after pancaking effect on halo velocities arises, it is gradually attenuated by 
gravitational force among halos, especially among very close halos. 
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THEORETICAL ANALYSIS 



In order to investigate the tendency of the anisotropic distribution of h alo velocity 
theore tically, we construct an analytic model using the simila.r met hodology to iLee fc Pen 
( 120011 ) ■ We begin on the Zel'dovich approximation (jZerdovichl Il970l ) 



X = q - L)(t)V^(q) 



'1^ 



which intrinsically estimates the formation of pancakes. Here, x is the Eulerian coordinate, 
q is the Lagrangian coordinate, \E'(q) is the linear velocity potential and D(t) is the growth 
factor, 



D = -n^H(a) 



da' 



[a'H{a')f' 



(2) 



By taking the time derivative of Zel'dovich approximation, we obtain equation of the velocity 
of cosmic particle smoothed on pancake scale. 



(3) 



Equation ([3]) is appropriate to investigate pancaking effect on halo velocity at early stage 
since equation ([T]) describes initial state of non-linear evolution. Also, we assume the density 
field is smoothed on the pancake scale, which is valid in this approximation because this 
approximation is available in the regime preceding the moment of the collapse to a pancake. 



Using a similar logic given by iLee fc PenI (120011 ) who found an expression for the align- 
ments between the tidal field and the halo position vectors, we take into account nonlinear 
effects into equation ([3]), and write an expression for the velocity alignments with the tidal 
field as 



5ij + r]K, 



= D'TikTki + 2DDT,kTki + 'TikTf 



- kj 



(4) 
(5) 



where T and T are tidal shear tensor and its time derivative, respectively. D is time 
derivative of D, which is 



where H'{a) 



D = -Qr 



dH{a) 



H'(a)d 



da' 



+ 



aH{a) 



(6) 



da 



{a'H{a')y am{a)\ ' 



Notice that r/ is a value in [0, 1]. If is 1, it represents a perfect correlation between local 
pancaking and halo velocity while if is 0, it means there is no correlation between them. 
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In velocity-velocity correlation, it is unlikely that the value of rj is close to unity since the 
velocity of a halo is not a conserved quantity and initial memory of local pancaking effect on 
velocity have been reduced. In order to show that the velocities of halos are not completely 
random but has some degree of alignments with the normal vectors to the local pancakes 
though The value of t] is expected to be very small, we have calculated the error on the value 
of T), that is, the standard deviation of t] for the case of no alignment. The forn iula for the 
error of r], a^^ is given as {A/lhNtotY^"^ where Ntot is the total number of halos ( Lee fc Pen 



2OOII ). With defined error, if 77 is larger than Sa^, one can say that there is a signal of true 



local pancaking effect. 

W e diagonalize equation 
fll996h 



by applying the relation referred by lBond. Kofman. &: Pogosyan 



where Aj and A, are an eigenvalue of the T and its time derivative, respectively. Thus, 
equation ^ is changed to 



The eigenvalues Aj are ordered by Ai > A2 > A3. 



212 



(9) 



To obtain the probability distrib ution, we assume probability density distribution of the 
alignment angles, cos 6*, is Gaussian (lLeell2004l ) 



p{v\'T) 



v/(27r)3det(M) io 



exp 



V dv, 



(10) 



where M, 



iVj\T). Finally, we can express the probability distribution of cos 6', p(cos6') 



IS, 



2 = 1 



2n 



sin^ 9 cos^ 
1 + 1]- 3t]Xki 



sin 6'sin^ 



cos^ 6 



Here, XKi = Ku/EiKl 



1 + r]- ?>7]Xk2 1 + - 3?7A/o 

(11) 



To compare with the results from simulation analysis, we use ACDM cosmological pa- 
rameters, h = 0.73, CTg = 0.9, Tcmbfl = 2.725K, flh = 0.045, = 0.25, which are the same 
values as the parameters used in Millennium Simulation. First, we calculate a by using 
following equation: 

^P{k)T\k)W^^{kR), (12) 



a\R,z) = D{z) 
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where D{z) is a grow th function, P ( k) is primordial power spectrum that is P{k) oc k, T{k) 
is transfer function (IBardeen et al.l Il986l ) and Wth is spherical top-hat window function. 
When evaluating equation (IT^ . we use Lp obtained by simulation analysis for R. Also, the 
redshift at which Lp is determined in the simulation is chosen as z value in equation f[T^ . 

Then we obtain Xki in equation (ITOl) . We decide Xki as 1 because Ai = 1 means the 
formation of the structure like a pancake in the Zel'dovich approximation, 



P 



P 



(1-Ai)(l-A2)(l-A3)^ 



(13) 



where p is the mean density of the universe and Ai, A2, A3 (Ai > A2 > A3) are the eigenvalues 
of the local tidal field, T. For \k2 and \k2.) we choose the most probable values of the 
following probability distribution. 



p(A2,A3|Ai = 1) =p(Ai = l,A2,A3)MAi = 1). 



(14) 



Her e, pfAi , A?, A3) and p(Ai ) are the probability distribution derived by iDoroshkevichI (119701 ) 
and iLee fc ShandarinI (119981 ) : 



p(Ai, A2, A3) 



3375 



8^5: 



p(Ai 



Una 
— \/2ti exp 

+ 3\/37rexp 



20— exp 
a 



■(Ai - A2)(A2 - A3)(Ai - A3) exp 



r2 



a' 



15/2 
2^ 



(15) 



9A| 
'2(t2 



5A| 
2a\ 

_15A| 
4(t2 



X2 



1 - 20^ +3V37rexp 



4a2 / I 2(T 



(16) 



We put these eigenvalues in equation ([9]) to obtain the normalized eigenvalues, Xki- 



Finally, we fit theoretical estimation of p(cos^), equation (ITT]) , to the numerical results 
as adjusting a correlation parameter r] in equation (fTTl) . We determine rj minimizing 
distribution 

(17) 



X 



where Xg^i and Xt^i are the values from simulation analysis and from theory, respectively and 
as,i is the Poissonian error from simulation analysis. The comparison numerical results with 
theory is shown in Fig. HI The filled circle points represent the numerically measured p{cos6) 
with Poissonian errors and solid line represents the analytic distribution. As can be seen in 
the figure, the simulation data points are in good agreement with our analytic formulation. 
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Also, as shown in Table [2l the strength of local pancaking effect on the halo velocity depends 
on the redshift. At 2; = 2, though halos are closely located at each other like Rc = 3h~^Mpc, 
7] is 0.011, which turns out to be 15.3(7,, (see Tabled]). Thus, even though the value of 77 is 
much less than unity, it is definitely a strong signal of velocity-pancake alignments, exceeding 
15 times the standard deviation, cr^. This strong signal which approximately amounts to 
15(7^ shows the existence of the velocity-pancake alignment. Even when i?c = at 2; = 2, 
is 0.003 which is 4cr^ which indicates that a pancaking effect appears even at Rc = 0. Then 
going to present epoch, the effect is reduced by gravitational attraction. Finally, at 2; = 0, 
when Rc is 3h~^Mpc, rj is almost and is also 0. In addition, as one can see the tendency 
in Fig. Hland in Table [2l at 2 = 2, pancaking effect becomes stronger as Rc is approximately 
5 — 6/i~^Mpc and then rj and gradually decreases when Rc is larger than 6h~^Mpc. In 
other words, the value of 1] and cr^ at 2; = 2 is the largest with i?c — 6. A similar tendency 
also appears at different redshifts, instead the range of Rc is shifted to larger scale caused 
by the expansion of the universe. Thus, we analytically find that there is a specific range of 
pancake scale where pancaking significantly affects halo velocity. In addition, we show how 
the pancaking effect changes with different redshifts. 

To find the tendency of the redshift dependence of 77, we fit 77 as a function of redshift 
to quadratic equation. The coefficients are shown in Table [31 Figure [5] plots the quadratic 
fitting function (dashed line) and compares it with the simulation data points (filled circles). 
As Fig. [5] shows, we note that rj is a. function of redshift, in a good agreement with the 
quadratic equation. Moreover, 77 is a function of the characteristic pancake scale, Lp. We 
also fit r]{Lp) to quadratic equation which can be seen in Fig. O The dashed line and the 
filled circles represent the fitting function and the data points, respectively. We can see a 
good match between the quadratic functional form and the simulation data. The best-fit 
coefficients are shown in Table |H 



4. DISCUSSION AND CONCLUSION 



We analyze the evolution of the pancaking effect by using millennium simulation. We 
obtain the signal of the strong alignment between the halo velocity and the normal to the 
local pancake plane at 2; ^ 2 when characteristic pancake scale is approximately 3 h~^Mpc. 
The alignment signal implies the first collapse of the protocloud may have occurred around 
2; ~ 2 and the size of the pan cake may have been around 3 h~^Mpc. This result is in 
agreement with iMo et al.l (120051 ) which reports pancakes form at 2; ~ 2 and their masses is 
approximately 5 x lO^^M©. 



Then as shown in Fig. [T] and Fig. [31 the alignment signal evolves. Pancaking effect 
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is weakened as going to low redshift, which means the pancaking effect on halo velocity is 
gradually attenuated by the gravitational attraction among the neighbor halos. In addition, 
pancake scale at low redshift is larger than that at high redshift. 

We expect the local pancaking effect could be detected in observation. It may be possible 
to achieve the signal of the velocity align ment using the velocity data reconstructed from 
2MASS redshift survey (jPirin et al.ll2006l ). We intend to measure the signal in observation 
for our future project. 

It is interesting to note that the pancakes are expected to have for med coincidentally 
at z when the star formation rate in massive galaxies is the highest (j Juneau et all 12006 



Feulner et al.ll2005l ). which is one of the phenomena referred to 'downsizing'. In 'downsiz- 
ing' scenario, the stars in more massive galaxies formed at higher redshift and those in less 
mass ive galaxies recently f ormed, where 'antih i erarch y' can be suggested. Recently, how- 



ever, 



Neistein et al.l (120061 ): iMouri fc Taniguchil (120061 ) show this phenomenon is natural in 



hierarchical clustering process. Nevertheless, a coincidence between the pancake formation 
epoch and the epoch of the heavily star-forming rnay irn ply the structure formation is not 
simply hierarchical but complicated (jCimatti et al.ll2006l ). 



In conclusion, we obtain the result that is the evolution of the velocity anisotropy 
induced by local pancake formation. This result would help to understand the unknowns in 
the galaxy formation. 



The Millennium Run simulation used in this paper was carried out by the Virgo Su- 
percomputing Consortium at the Computing Centre of the Max-Planck Society in Garching 
(http: //www.map-garching.mpg.de/millennium|). We are thankful to L. Gerard for the halo 
catalogs of the Millennium Run simulation and to V. Springel for useful discussion. We also 
appreciate the anonymous referee for helpful suggestions. This work is supported by the 
research grant No. ROl-2005-000-10610-0 from the Basic Research Program of the Korea 
Science and Engineering Foundation. 



- 10 - 
REFERENCES 

Aragon-Calv, M. A. ,van de Weygaert, R. Jones, B. J. T., & van der Hulst, J. M. 2007, ApJ, 
in press ( astro-ph/0610249) 

Bardeen, J.M., Bond, J.R., Kaiser, N., & Szalay, A.S. 1986, ApJ, 304, 15 

Bond, J., R., Kofman, L., & Pogosyan, D. 1996, Nature, 380, 603 

Brunino, R., Trujillo, I., Pearce, F. R., & Thomas, P. A. 2007, MNRAS, in press 
( [astro^/0609629| ) 

Cimatti, A., Daddi, E., & Renzini, A. 2006, A&A, 453, 29 

Croton, D. J., Springel, V., White, S. D. M., De Lucia, G., Frenk, C. S., Gao, L., Jenkins, 
A., Kauffmann, G., Navarro, J. F., & Yoshida, N. 2006, MNRAS, 367, 864 

Doroshkevich, A. G. 1970, Astrofizika, 3, 175 

El-Ad, H., & Piran, T. 1997, ApJ, 491, 421 

Feulner, G., Gabasch, A., Salvato, M., Drory, N., Hopp, U., & Bender, R. 2005, ApJ, 633, 
L9 

Glazebrook, K. et ah 2005, Nature, 430, 181 

Hahn, O., Porciani, C., Carollo, M. C., & Dekel A. 2007, MNRAS, 375, 489 

Hoyle, F., & Vogeley, M. S. 2002, ApJ, 566, 641 

Juneau, S. et ah 2005, ApJ, 619L, 135 

Lee, J. & Shandarin, S. F. 1998, ApJ, 500, 14L 

Lee, J., & Pen, U. L. 2001, ApJ, 555, 106 

Lee, J., & Pen, U. L. 2002, ApJ, 567, Llll 

Lee, J. 2004, ApJ, 614, LI 

Mo, H. J., Yang, X., van den Bosch, F. C., & Katz, N. 2005, MNRAS, 363, 1155 
Mouri, H., & Taniguchi, Y. 2006 A&A, 459, 371 
Noh, Y., & Lee, J. 2006, ApJ, 652, L71 



- 11 - 



Neistein, E., van den Bosch, F. C, & Dekel, A. 2006, ApJ, 372, 933 

Erdogdu, P., Lahav, O., Huchra, J. P., CoUess, M., Cutri, R. M., Falco, E., George, T., 
Jarrett, T., Jones, D. H., Macri, L. M., Mader, J., Martimbeau, N., Pahre, M. A., 
Parker, Q. A., Rassat, A., & Saunders, W. 2006, MNRAS, 373, 45 

Shandarin, S. F., Melott, A. L., McDavitt, K., Pauls, J. L., & Tinker, J. 1995, Phys. Rev. 
Lett., 75, 7 

Spergel, D. N. et al. 2003, ApJ, 148, 175 

Springel, V. et al. 2005, Nature, 435, 629 

Zel'dovich, Ya. B. 1970, A&A, 5, 84 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 12 - 



1 .04 


z = 0.0_ 


z=1.0_ 


z = 2.0_ 


1 .02 




1 1 \ I 


1 I 1 


1 .00 


- - 


1--- 


--\-v-i-\-\ 1-^-1-^- - 


^ i 




1 M 1 t t 1 1 


0.98 


V 


I 


- 


0.96 






- 


1 .04 


- 




1 


1 .02 




1 ill" 


1 

} 1 r 


1 .00 


I 1 , 1 

4 iH 1 1 , ,- 


--i-i - 


- — 


1 ' t 1 • 




n 1 1 ' 


0.98 


J I 




, 1 1 , 1 


0.96 






- 


1 .04 








' .02 








' .00 














0.98 








0.96 









0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 

COS0 COS0 COS0 



Fig. 1. — Probability density distributions of the cosines of the relative angles between the 
field halo velocities and the directions normal to the local planes enclosing the two neighbor 
field halos for the cases of the distance threshold: Rc = 0,3,6 /i~^Mpc at different redshift: 
z = 0, 1, 2 from Millennium Run simulation. The errors are Poissonian and the horizontal 
dotted lines represent no alignment. Note that halos which contain more than 50 particles 
are determined 
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Fig. 2. — The same figure as Fig. [T]but using the halos which contain more than 30 particles. 
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Fig. 3. — The average of the cosines angles (cos^) between the halo velocities and the 
local pancake planes whose definition is the same as Fig. [1], corresponding to the change of 
characteristic linear size of the pancake, Lp at four different redshift: ^ = 0, 0.5, 1, 2. The 
horizontal dotted lines represent no correlation. 
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Fig. 4. — Probability density distributions of the cos ^s. Filed circle with Poissonian errors 
represent the results from the simulation while solid line stand for the analytic prediction 
with correlation parameter r\. The horizontal dotted lines correspond to no correlation. 
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Fig. 5. — r/ as a function of redshifts. Filed circle represents the numerical measurement in 
Millennium Simulation. Dashed line is a quadratic fitting function at different R^. scale. 
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Fig. 6. — The same as Fig. [5] but as a function of characteristic pancake scale, Lp. 
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Table 1. left: Total number of identified halos whose contain more than 50 particles in 
each redshift catalog, right: The number of field halos 



z 


N total 


Nfield I 


'Infield 





6126689 


502108 


(8.2) 


0.5 


6631145 


542554 


(8.2) 


1 


6817225 


558138 


(8.2) 


2 


6222187 


516774 


(8.3) 
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Table 2. Correlation parameter rj corresponding to Rc and redshift, z with its standard 
deviation for the case of no alignment. The formula for the error of V, o'rj is given as 
{4:/15Ntoty^'^ where Ntot is the total number of halos (ILee fc Penll200l[ ). 



Rc 






z 







0.5 


1 


2 


3 


0.000(0.0(7,,) 


0.003(4.3a^) 


0.004(5.8(7^) 


0.011(15.3(7^) 


6 


0.005(6.9(7,,) 


0.008(11.4a^) 


0.009(13.0(7^) 


0.013(18.1(7,,) 


9 


0.007(9.6cr^,,) 


0.009(12.8(7,,) 


0.011(15.9(7^) 


0.012(16.7(7,,) 
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Table 3. Fitting parameters for 77 as a function of redshift, z. Fitting functions are 
quadratic function, r]{z) — az^ + 62; + c, (10"^) 



z 


a 


h 


c 





1.36 


-1.28 


0.08 


0.5 


1.00 


3.30 


0.30 


1 


-0.36 


4.58 


5.22 


2 


-1.36 


5.28 


6.92 
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Table 4. The same as Table [3] but as a function of linear characteristic pancake scale, Lp, 







+ bLp + c, 


(10-3) 


z 


a 


b 


c 





-0.31 


4.73 


-9.71 


0.5 


-0.40 


5.36 


-9.26 


1 


-0.37 


4.76 


-6.62 


2 


-0.30 


3.39 


2.58 



